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An Efimov trimer near the atom-dimer threshold can increase the atom loss rate in ultracold 
trapped atoms through the avalanche mechanism proposed by Zaccanti et al. A 3-body recom- 
bination event creates an energetic atom and dimer, whose subsequent elastic collisions produce 
additional atoms with sufficient energy to escape from the trapping potential. We use Monte Carlo 
methods to calculate the average number of atoms lost and the average heat generated by recom- 
bination events in both a Bose-Einstein condensate and a thermal gas. We take into account the 
energy-dependence of the cross sections and the spatial structure of the atom cloud. We confirm 
that the number of atoms lost can be much larger than the naive value 3 if there is an Efimov trimer 
near the atom-dimer threshold. This does not produce a narrow loss feature, but it can significantly 
affect the determination of Efimov parameters. 
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Particles with short-range interactions and an S-wave 
scattering length a that is large compared to the range 
have universal low-energy properties that depend on a 
but not on other details of the interactions [1] . Universal- 
ity provides connections between various fields of physics, 
including atomic and molecular, condensed matter, nu- 
clear, and particle physics. This has stimulated theoret- 
ical progress in few-body physics. Further stimulus has 
come from the use of ultracold trapped atoms to study 
reaction rates at extremely low energies. 

Since particles with large scattering length are essen- 
tially indivisible at low energies, we will refer to them 
as atoms. In the 2-atom sector, the universal proper- 
ties for a > include the existence of a shallow dimer 
(weakly-bound diatomic molecule). In many cases, in- 
cluding identical bosons, the universal properties in the 
3-atom sector include the existence of a sequence of uni- 
versal triatomic molecules called Efimov trimers [2J. In 
the zero-range limit, the spectrum of Efimov trimers is in- 
variant under discrete scale transformations [3 . For iden- 
tical bosons, the discrete scaling factor is approximately 
22.7. Reaction rates among three low-energy atoms also 
respect discrete scale invariance [3]. We refer to univer- 
sal few-body phenomena with discrete scaling behavior 
as Efimov physics. 

Ultracold trapped atoms provide an ideal laboratory 
for studying Efimov physics, because a can be controlled 
experimentally using Feshbach resonances. The simplest 
probes of Efimov physics are loss features: local max- 
ima and minima in the atom loss rate as functions of a. 
The most dramatic signature of an Efimov trimer is the 
resonant enhancement of the 3-body recombination rate 
when there is an Efimov trimer near the 3-atom thresh- 
old [5]. The first observation of such a loss feature in 
an ultracold gas of 133 Cs atoms [6] revealed a line shape 
consistent with universal predictions [7]. 



In a mixture of atoms and shallow dimers, a narrow 
loss feature can also be caused by an Efimov trimer near 
the atom-dimer threshold. We will refer to a scattering 
length a* for which an Efimov trimer is exactly at the 
threshold as an atom-dimer resonance. For a near a*, 
there is resonant enhancement of both the elastic scat- 
tering of an atom and the shallow dimer and their inelas- 
tic scattering into an atom and a deep dimer (strongly- 
bound molecule). The large binding energy of the deep 
dimer gives the outgoing atom and dimer large enough 
kinetic energies to escape from the trapping potential. 
The resulting peak in the atom loss rate near a* was first 
observed in a mixture of 133 Cs atoms and dimers [5]. 

There have also been observations of enhanced loss 
rates near a* in systems consisting of atoms only. Zac- 
canti et al. observed a narrow loss peak near the predicted 
position of an atom-dimer resonance in a Bose-Einstein 
condensate (BEC) of 39 K atoms [S]- They also observed a 
loss peak in a thermal gas near the next atom-dimer res- 
onance, at a scattering length larger by a factor of about 
22.7. Pollack et al. observed a loss peak near the pre- 
dicted position of an atom-dimer resonance in a BEC of 
7 Li atoms [TU]. Machtey et al. observed such a loss peak 
in a thermal gas of 7 Li atoms [IT]. These loss features 
near the atom-dimer resonance are puzzling, because the 
equilibrium population of shallow dimers is expected to 
be negligible in these systems. 

Zaccanti et al. proposed an avalanche mechanism for 
the enhancement of the atom loss rate near a* in systems 
consisting of atoms only [9] . The loss of atoms is initiated 
by a 3-body recombination event that produces an atom 
and a shallow dimer with kinetic energies large enough to 
escape from the trap. If they both escape, there would 
be 3 atoms lost. If the dimer instead scatters inelasti- 
cally, the scattered atom is also lost. However the dimer 
can undergo multiple elastic collisions before ultimately 
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FIG. 1: (Color online) Cross sections for elastic (upper panel) 
and inelastic (lower panel) atom-dimer collisions for 77, = 0.2. 
The straight lines are elastic cross sections for atom-atom 
collisions. The three curves for each case (in order of generally 
increasing cross sections) are for the first scattering, a typical 
second scattering, and after many elastic scatterings. The 
vertical dotted lines mark the positions of a* and 22.7 a*. 



escaping or suffering an inelastic collision, and it may de- 
liver enough energy to the scattered atoms to allow them 
to escape from the trap. These atoms may also undergo 
multiple elastic collisions, resulting in an avalanche of 
additional lost atoms. Near a*, the resonant enhance- 
ment of the atom-dimer cross sections increases both the 
probability for the dimer to initiate an avalanche and 
the probability for an inelastic collision. The resulting 
increase in the number of atoms lost per recombination 
event could produce an observable loss feature. 

In this paper, we present a quantitative analysis of 
the avalanche mechanism for atom loss. We use Monte 
Carlo methods to generate avalanches of atoms initiated 
by recombination events with the appropriate probability 
distribution. We calculate the average number of atoms 
lost and the average energy converted into heat from an 
avalanche. We use the results to calculate the atom loss 
rate constant for both a Bose-Einstein condensate and a 
thermal gas of trapped atoms. 

We consider identical bosons of mass m with a large 
positive scattering length a. The 2-body physics in our 
Monte Carlo model for the avalanche mechanism consists 
of the binding energy = h 2 /ma 2 for the shallow dimer 
and the cross section <taa = 87ra 2 / (1 + a 2 k 2 m ) for elastic 
atom-atom scattering with center-of-mass wavenumber 



k cm . The 3-body physics in our model consists of the 
rate constants a s haiiow and cvdeep for 3-body recombina- 
tion at threshold into the shallow dimer and into deep 
dimers and the cross sections and a^p for elastic 
and inelastic atom-dimer scattering. In the zero-range 
limit, these reaction rates are determined by a and two 
Ehmov parameters: the atom-dimer resonance a* and a 
dimensionless parameter 77* that controls the decay width 
of an Efimov trimer [7] . Analytic expressions for a s haiiow 
and ttdeep are given in Ref. [1 . They are more conve- 
niently expressed in terms of a parameter a*o that differs 
from a* by a universal ratio: a*o/a* = 4.4724. The 5 
digits of accuracy in the ratio are obtained by combining 
universal results from Refs. [T21|T3] . Parameterizations of 
<j^d and cr^ 1 ^ from the atom-dimer threshold up to the 
dimer-breakup threshold can be obtained from Ref. [T]. 

We make several simplifying approximations in our 
model. We ignore the effects of potential energies on the 
cross sections for the dimer and the atoms. In the case 
of a BEC, we also ignore mean-field energies. We also 
approximate the trajectories of the atoms and the dimer 
between collisions as straight lines. In the first collision 
of the atom or dimer from the recombination event with 
a stationary atom, the center-of-mass wavenumber fc cm is 
l/(y3a) or 2/(3\/3a), respectively. In the second colli- 
sion with a stationary atom, the typical /c cm is smaller by 
a factor of 1/ \/2 for the atom and V5/3 for the dimer. As 
the number of elastic collisions increases, fc cm decreases 
towards 0. The universal cross sections for the first col- 
lision, a typical second collision, and after many elastic 
collisions (k cm — > 0) are shown in Fig. [I] for 77* = 0.2. 
At k cm = 0, 

a AD an d ^cm^AD nave dramatic peaks with 
maxima near a* and 22.7 a*. The elastic cross section 
also has deep minima near 0.38 a* and 8.6 a*. For the 
first few collisions, ct^d still has deep minima but there 
are no dramatic peaks. 

The experimental inputs in our Monte Carlo model 
are the number Nq of trapped atoms, the frequencies v Xl 
u y , and v z of the harmonic trapping potential, the tem- 
perature T of the atoms, the scattering length a (which 
can be controlled by varying the magnetic field near a 
Feshbach resonance), and the trap depth Strap: which 
should be much larger than the energies of the trapped 
atoms. Atoms and dimers that reach the edge of the 
atom cloud are assumed to be lost if their energies ex- 
ceed -Etrap and 2i5 t rap, respectively. The role of the re- 
maining experimental inputs is to determine the number 
density n(x, y, z) of the trapped atoms. We consider two 
simple cases: a BEC of atoms at zero temperature in the 
Thomas-Fermi limit and a thermal gas of atoms above 
the critical temperature for BEC. The rate at which the 
number N of atoms in a thermal gas decreases due to 
3-body recombination is 

dN/dt = - ((Most) "shallow + 3a d co P ) (n 2 )N, (1) 
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where (n 2 ) and (Most) are spatial averages weighted 
by n(r) and n 3 (r), respectively. The right side must 
be multiplied by 1/3! if the system is a BEC. Atoms 
produced by the avalanche that have energy less than 
Strap can never escape from the trapping potential and 
their kinetic energy will ultimately be transformed into 
heat. The recombination heating rate in a thermal gas is 
(^hcat)a s haiiow(« 2 )^V, where (S hcat ) is the average heat 
from a single avalanche. 

The development of an avalanche can be decomposed 
into discrete steps corresponding to the recombination 
event and the subsequent scattering events. For each 
event, the subsequent state of the avalanche has a sim- 
ple probability distribution, (a) The position (x, y, z) 
of the recombination point has a distribution propor- 
tional to n 3 (x,y,z). (b) The momenta of the outgo- 
ing particles from an event have a distribution that is 
isotropic in the center-of-mass frame, (c) An atom or 
dimer flies beyond the edge of the atom cloud with prob- 
ability exp(— a fnd£), where a is a a a for an atom and 

a AD + a ad f° r a dimer and where Jn d£ is the column 
density integrated from the position of the previous colli- 
sion out to infinity along a straight path in the direction 
of the momentum. If the random number determines 
that an atom or dimer fails to reach the edge of the atom 
cloud, it is also used to determine the position where it 
scatters, (d) Given that a dimer scatters, it scatters in- 

elastically with probability o'ad/^Ad +°ad)- All these 
simple probability distributions together determine the 
probability distribution of avalanches. 

We generate avalanches with the appropriate proba- 
bility distribution using a Monte Carlo method that pro- 
duces a binary tree whose nodes represent events. The 
branches represent the two outgoing particles from each 
event. There are also terminal nodes that correspond 
to atoms and dimers whose ultimate fate is determined. 
The conditions for a terminal node depend on the kinetic 
energy S of the particle, (a) If an atom has E < St rap , it 
remains trapped, (b) If an atom that reaches the edge of 
the atom cloud has E > Strap, it is lost- (c) If a dimer has 
an inelastic collision, both it and the scattered atom are 
lost, (d) If a dimer that reaches the edge of the cloud has 
E > 2S trap , it is lost, (e) If a dimer that reaches the edge 
of the cloud has E < 2S trap , it will return to the cloud 
and eventually suffer an inelastic collision. The terminal 
nodes give contributions to the number of atoms lost and 
to the heat of the remaining atoms. Adding these con- 
tributions from all the terminal nodes, we get Most and 
Ehc&t for the avalanche. We calculate (Most) and (Sh ea t) 
by averaging over many avalanches. More than 100,000 
avalanches are sometimes required to get smooth results 
for (Most) and (Sheat) as functions of a. 

Zaccanti et al. developed a simple probabilistic model 
for the avalanche process [9|. In the Zaccanti model, the 
avalanche is reduced to a discrete sequence of dimer scat- 
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FIG. 2: (Color online) Rate constant L3 for 3-body recom- 
bination in 7 Li atoms. The data are for the thermal gas 
from Ref. |11| . The lower solid curve is the naive result for 
(Most) = 3 with a* = 282 a and r\, — 0.2. The upper solid 
(red) and dashed (blue) curves are predictions of our Monte 
Carlo model for the thermal gas in Ref. [TT] and for the BEC 
in Ref. [TO]. The vertical dotted line marks the position of a*. 

tering events. A variable number of elastic collisions is 
followed either by the escape of the dimer from the trap 
or by a final inelastic collision. There is one lost atom for 
each elastic collision up to a maximum number that is 
determined by the trap depth S trap . The relative proba- 
bility for each sequence of scattering events is determined 

by the mean column density and by c^d ana - ^cm.o'AD- 
The Zaccanti model is greatly simplified in several ways 
compared to our Monte Carlo model: (a) The spatial 
structure of the avalanche is ignored, (b) Elastic scatter- 
ing of the atoms is not considered, (c) The energy de- 
pendence of Cad ano - tmfjo is not taken into account. 
Zaccanti et al. used their model to calculate (Most) for 
their experiment with 39 K atoms [5]. It predicts that 
(Most) increases from its background value of 3 to about 
13 near the atom-dimer resonance. The resulting predic- 
tion for the atom loss rate agrees qualitatively with the 
loss feature they observed near 30.4 ao. 

Machtey et al. developed an alternative probabilistic 
model for the avalanche process [2]. They made the 
same simplifications as in the Zaccanti model, but they 
used different probabilities for the sequences of scatter- 
ing events. Machtey et al. did not introduce the trap 
depth Strap, so they could not calculate (Most)- Instead 
they calculated the average number M of dimer collisions. 
They suggested that the loss feature near a* might be as- 
sociated with the maximum of N. 

In the experiment with 39 K atoms in Ref. jS], the loss 
feature near 30.4 ao is at a scattering length that may 
be too small for universal predictions to be reliable. We 
therefore focus on the experiments with 7 Li atoms. In 
Ref. [10], the data near the atom-dimer resonance were 
obtained using a BEC of 7 Li atoms in the |1,+1) by- 
perfine state with Strap ~ 0.5 //K. They measured the 
rate constant L3 = 3(a sria iiow + c*deep) as a function of 
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FIG. 3: (Color online) Average number (Most) of atoms 
lost (upper panel) and average heat (Ehc&t) generated (lower 
panel) by an avalanche from 3-body recombination in 7 Li 
atoms. The dashed (blue) and solid (red) curves are the pre- 
dictions of our Monte Carlo model with a* = 225 ao and 
?7» = 0.2 (lower curve) or rj, — 0.04 (higher curve) for the 
BEC in Ref. [TO] and for the thermal gas in Ref. [EE]. The 
vertical dotted lines mark the positions of o* and 22.7 a*. 



a. The best fit to L3 near an interference minimum gave 
the Efimov parameters a*o = 2672 ao and 77* = 0.039. 
The accuracy of the determination of a as a function of 
the magnetic field has been improved and a reanalysis 
of that data is underway |15j . The parameter a^o will 
change significantly, but 77* will not. In Ref. |16j . L3 was 
measured as a function of a using a thermal gas of 7 Li 
atoms in the |1,+1) hyperfine state with T 1.4 /iK 
and Strap ~ 7 /xK. The data are shown in Fig. [2j The 
best fit to L3 gives a*o = 1260 ao and 77* = 0.188. Using 
the universal ratio for a*o/a*, we obtain the prediction 
a* = 282 ao. The curve for a* = 282 ao and 77* = 0.2 
shown in Fig.[2]provides a good fit over most of the range. 
The data in Fig. [2] for a below 220 ao were presented in 
Ref. [TT], revealing the narrow loss feature near 200 ao- 

In Fig. [3j we show the predictions of our Monte Carlo 
model for (Most) and (Sheat) as functions of a for both 
7 Li experiments. We set a* — 282 ao and we consider 
both 77* = 0.2 and 77* = 0.04. There are discontinuities 
in (Most) and (-Ejieat) when the energies of the atom and 
the dimer from the recombination event are equal to Strap 
and 2£;trap, respectively. Both (Most) and (E hcat ) have 
broad peaks whose maxima occur well above a* . This is 
because the first few elastic collisions of the dimer from 



the recombination event are the most important for gen- 
erating an avalanche. As shown in Fig. [TJ the maxima in 
the elastic cross sections for the first few collisions occur 
well above a*. Decreasing 77* by a factor of 5 does not 
give dramatic changes in (Most) and (-Eheat)- 

If the avalanche mechanism is taken into account, as in 
Eq. 0, the rate constant is L 3 = (Most) "shallow + 3a dccp . 
In Fig. [2] we show the predictions for both experiments 
for L 3 as a function of a using a* = 282 ao and 77* = 0.2. 
If 77* were decreased to 0.04, the changes in the curves 
would be obvious only near a*o ~ 1260 ao, where the 
minimum would be about a factor of 5 below the data. 
There is no narrow loss feature near a* , but instead there 
is a broad enhancement in L3 in the region between a* 
and a*o- The enhancement is large enough that it could 
affect the fitted values of the Efimov parameters. We are 
unable to get a narrow loss feature like that in the data 
in Fig. [2] for any values of the parameters. 

Our results for the avalanche mechanism also suggest 
that the effects of heating should be taken into account 
differently in the analysis of atom loss data at positive 
scattering lengths. In addition to the recombination 
heating associated with -E^eat, one must take into ac- 
count the disappearance of the three atoms that undergo 
recombination |17j . In a thermal gas, the resulting rate 
of change in the temperature is 

dT/dt = ( [<£ heat )/(3fcT) + (Most) - 2] 

^shallow 

+a dccp )(n 2 )T. (2) 

This has the same dependence on T as that assumed 
in Ref. [TT]. In the coefficient of cv s haiiow, the 1/kT 
term is determined by (E\ ost ) and the additive constant 
(Most) — 2 can differ from the naive value 1. The effects 
of heating are usually taken into account by using the 
coupled equations for dN/dt and dT/dt to extrapolate to 
the initial value of dN/dt [TT]. The (Most) and (Shcat) 
terms in Eqs. ([!]) and ([2| may have a significant effect on 
this extrapolation. 

We have found that the avalanche mechanism does not 
produce a narrow loss feature near the atom-dimer reso- 
nance, but instead a broad enhancement of L3 between 
a* and a*o- This can have a significant effect on the 
determination of Efimov parameters from data at posi- 
tive scattering lengths. The heating from the avalanche 
mechanism can also be important in the extrapolation 
to the initial atom loss rate. All experiments on Efimov 
loss features at positive a should probably be reanalyzed 
to take into account the effects of the avalanche mecha- 
nism. If the observed loss features near a* in 39 K and 
7 Li atoms survive such a reanalysis, the mechanism for 
these loss features will remain a puzzle. 
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